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In this work, the integrability of the one dimensional chiral Hubbard model is 



CN . discussed in the limit of strong interaction U = oo. The system is shown to be 

> 



integrable in the sense of existence of infinite number of constants of motion. The 
system is related to a chiral Kondo lattice model at strong interaction J = +cx3. 
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Hubbard model has been of considerable interests due to its possible relevance to the high 
temperature superconductivity. The one dimensional Hubbard model with nearest neigh- 
boring hopping was solved by Lieb and Wu in 1968 with the help of Bethe-ansatz [|l]-0]. It 
was shown that the system exhibits a metal-insulator phase transition at half-filling even 
for arbitrarily small interaction. Away from half-filling, the low lying excitations of the 
system have been classified as Luttinger liquid like in the sense of Haldane. About two 
years ago, a one dimensional SU{2) Hubbard model with only relativistic right movers was 
introduced, which reduces, at half-filling and large but finite on-site energy, to the SU{2) 
Haldane- Shastry spin system with 1/r^ exchange interaction [Q. Using finite size diagonal- 
ization result and the information provided by some special cases, an effective Hamiltonian 
was proposed, which was used to provide the full energy spectrum and the thermodynamics 
for any on-site energy. With the help of the effective Hamiltonian, it was found that at 
T = and half-filling, there exists a critical value Uc at which metal-insulator phase tran- 
sition occurs in the system ^. It was conjectured that the system is completely integrable 
for any on-site energy. However, a proof for the conjectured energy spectrum and the ther- 
modynamics, as well as the structure of the wavefunctions at finite t/, is still unknown. The 
integrability of this Hubbard model remains an open problem. 

In the strong interaction limit t/ = oo, it was discovered that the Gutzwiller-Jastrow 
product wavefunctions are eigenstates of the chiral Hubbard model, both in the SU{2) 
case and in the SU{N) case 0. In fact all eigenstates can be expressed in terms of more 
generalized Jastrow product wavefunctions. Furthermore, the SU{2) energy spectrum |^ is 
the same as conjectured by Gebhard and Ruckenstein in their original work 0]. In this work, 
we would like to discuss the integrability of the system in the limit where the interaction 
between the electrons is infinitely strong U = oo. Using simple argument, we shall exhibit 
an infinite number of constants of motion, showing that the system is integrable. 

In the second part of our paper, we consider a one dimensional chiral Kondo lattice. 
The conduction band has only right moving electrons, and the electrons interact with each 
localized impurity moment through exchange interaction. We identify the chiral Kondo 



lattice at J = +00 with the Hubbard model which has been previously studied. With this 
identification, the full energy spectrum, the wavef unctions, the thermodynamics and the 
integrability of the system can be obtained for the Kondo lattice in this limit. In particular, 
various correlation functions between the electrons and the impurity spins can be computed 
exactly for this system. 

The chiral Hubbard model Hamiltonian is defined on a one dimensional lattice of length 
L: 

N 

H= J2 Y^i^ij^tcja) + U J2 H ^i^^i<T', (1) 

l<ijtj<L(7=l l<i<La^a' 

where the hopping matrix element is given by t^n = (— «t)(— l)(™^"')[(L/7r) sin(7r(m. — 
n)/L)]^^. For this SU{N) system, the spin of the electrons can take values from 1 to 
A^. In the case of SU{2), it is the Hamiltonian introduced by Gebhard and Ruckenstein 
1^. In the strong interaction limit U = cxo, at each site there is at most one electron, the 
number of holes N^ and the number of the electrons on the lattice N^ are conserved quanti- 
ties. To rewrite the Hamiltonian in more convenient form, we perform the following unitary 
transformation T: 

J _. ('_1 ^3; -TTiz/L t 

under which the original Hamiltonian becomes 

H^H = 2tPa[ E Ht-^Ac.APgj. (3) 

with Pgi the Gutzwiller projector, making sure that there are no double or multiple occu- 
pancies, while Zx = e^'^"^^, with x = 1, 2, ■ ■ ■ , L. 

In the Hilbert space of no double or multiple occupancies, the electron fields can be 
rewritten with the superalgebra representations, 

PGi^)4.PGi^) = flk, 

PG{^)c^.PG{^) = h\h,, (4) 



where the / fields are fermions, the b fields are bosons, with the constrain J2a fLfirr+blbi = 1. 
Pg(0 is the Gutzwiller projector operator on the site i for the electron operators c and c^. 
Any state vector of the Hilbert space can be written as 

10 >= E 0(w},M)n/.unHjo>, (5) 

{xa},{y} i=l j=l 

where the amplitude is symmetric in the coordinates {y} of the b bosons, while anti- 
symmetric when exchanging the spin and positions Xiai,XjCrj of two / fermions. Here, 
illy '?2, ■ ■ ■ ; 1l) = (a^i, ■ ■ ■ , X]\[^,yi, ■ ■ ■ , yj\fj span the full chain. The eigenenergy equation of 
the system can then be written in the first quantized form as follows: 
9fiT ^e ^ z 1 

^E E Tv^y-^M^^-o E ^^(M,W) = ^0(W>W), (6) 

where Zi = ^'^'^'nl^ ^ the operator Mj^ exchanges the position variables g, and g^, the operator 
Pjj exchanges the / fermion spin variables Oi and Oy These two operators commute with 
each other, as they act on different groups of variables of the wavefunction. 

Following the ideas of Ref. |TT]], we define generalized momentum operators IIj, with 



z = l,2,---,iV„ 

n,= 5: \/,,M,„ (7) 

with Vij = Zj/{Zi — Zj). With this special form of V^jj, the generalized momentum operators 
satisfy the commutation relation 

[ni,n,] = Mi,ni-n,Mi,-, (8) 

where 1 < i 7^ j < iVg. We then introduce the following hermitian operators 

^n = En^ (9) 

s=l 

where n = 0, 1, 2, ■ ■ ■ , cxd, and the sum s is over the electrons, i.e., from 1 to N^.. In particular, 
the Hamiltonian is given by 

^ = ^[^1-^ E m- (10) 



Using the comniutation relations Eq. (P) it can be shown that all the operators A'^s commute 
with each other. Furthermore, the action of A^ on some amplitude (p does not change the 
symmetry properties, i.e., the resultant wavefunctions remains symmetric in the b boson 
positions, and antisymmetric when exchanging any pair of the / fermion positions and spins 
simultaneously. It is straightforward to prove the following relations 

M,,K0]= [A„0], iV, + 1 < z ^ J < L 

M,,/^;j[A„0]= (-i)[A„0], 1 < z ^ J < iv,. (11) 

Since all the operators A'^s commute with the Hamiltonian given by the Eq. ([T0|) , we thus 
have an infinite set of conserved physical quantities of the system, showing that the system 
is indeed completely integrable. With these A'^s, we can construct corresponding quan- 
tities written in second quantized language, which commute among themselves and with 
the Hamiltonian H. Carrying out the unitary transformation T~^, it is straightforward to 
convert them, so that the resultant quantities are constants of the original Hamiltonian H. 
One interesting observation is that the mutually commuting hermitian quantities A^s 
are also the invariants of the long range supersymmetric t-J model [Q, i.e., [An, Ht-j] = 0. 
As noted previously, the physical quantities. In = Y.i=i n", with n = 0, 1, 2, ■ ■ -, commute 
with each other, and with the supersymmetric t-J model Hamiltonian. The two fami- 
lies of quantities {/«}, {An} are independent of each other, in the sense that we can not 
write any member of one group in terms of a linear combination of members of the other 
group. Furthermore, by explicit computation, one can show that they do not commute, e.g., 
[y4„, Jm] 7^ 0. Therefore, it is clear that the previous family of conserved quantities {In}, 
although providing a proof of the integrability of the t-J model, does not exhaust all the 
constants of motion. The two infinite symmetries of the system do not commute with each 
other ip. It would be very interesting to find larger group of mutually commuting constants 
of motion. This is actually one of the most fundamental questions encountered in study of 
the quantum integrable systems. That is, when one has found an infinite number of simul- 
taneous constants of motion, showing that the system is integrable, it is still not sure that 



the infinite set would contain all possible simultaneous constants of motion of the system. 
In the following, we will show how the results of the chiral Hubbard model may be 
generalized to a one dimensional Kondo lattice model. The Kondo lattice model has been 



an interesting model for the study of heavy fermion systems |T0[ . In this model, the system 
has an array of localized impurity moments, the conduction electrons interact with the local 
moments through spin exchange. In general, the conduction band is best described by 
the tight banding picture, in which we have both right and left movers. However, in the 
following, we assume that the electrons propagate in only one direction. This chiral Kondo 
lattice model is defined on a one dimensional lattice, with the Hamiltonian: 

h=j:j: <k)cic,. + JEcL^Q, ■ Sf{^, (12) 

where the conduction band spectrum is e{k) = —tk, with k = ^^, — (L — l)/2 < K < 
{L — l)/2, in the momentum space. J is the coupling constant between the local impurity 
moments and the conducting electrons. The local moments are described by the spin 1/2 
operators, that is, [Sj{k), S^{k)] = iSf{k) (plus two other commutation relations obtained 
by the cyclic permutations of x,y,z), with the relation Sj{k) = 3/4, for all the sites k = 
1, 2, ■ ■ ■ , L. The Hamiltonian may also be written in the following way: 

H= Y. 12 ^iJ^aCj'^ + JY.^a^Cif^-Sfii), (13) 

i<i^i<i(7=Ta i=i 

where t^„ = (-it)(-l)(™-")[(L/7r) sin(7r(m - n)/L)]-i. 

When the interaction of the electron and the impurity is very strong, i.e. J = +00, we 
can map the system onto the above chiral Hubbard model with infinity repulsion. Indeed, 
when there are N^ electrons on the lattice L, with A^e ^ L, then each electron will attempt 
to form a singlet with the impurity spin at each site, to lower the energy of the system as 
much as possible, and some unpaired impurity spins are left over on the lattice. The Hilbert 
space at each site can be either a unpaired impurity spin or a singlet of electron-impurity 
bound state. Due to the hopping of the conduction electrons, the singlets can hop on the 
lattice. In this case, the basis vectors can be written as 



-N, 



11(1 -^7.^; 



.1=1 



c^ c^ ■ ■ -c^ 10 > 



\cri,<^2,---,(3l,--- ,P2,---,<^L-Ne >, (14) 



where the singlets are located at positions {x} = {xi < X2 < ■ ■ ■ < xn^), the unpaired 
impurity spins {cri,<y2,- ' ' y'^L-Ne) S'-re positioned at sites {y} = {yi < y2 < • • • < Vl-n^)- 
Here, the operator P^^p^ permutes the spin indices 7j and /9j, to form a singlet of electron 
and impurity at site Xj. With P the projector onto this subspace, the Hamiltonian takes 
the form: 

H = PHP = PTP + c, (15) 

where T is the kinetic energy of the conduction electrons, the infinite constant c = {—J/4)Ne 
only shifts the origin of the energy of the system, a reference energy which is unimportant 
physically. In the space where the z-component of the total spin is fixed, that is, 5*2 = M, 
the number of the unpaired up-spin impurities is A = M + (L — Ne)/2, the number of the 
unpaired down-spin impurities is -B = —M + (L — Ne)/2. C^" x C^_j^^ is the size of the 
Hilbert space. Any eigenstate of the Hamiltonian Hi = PTP can be written as a linear 
combination of the basis vectors, 

\(t)>=Y,C{a)\a> . (16) 

a 

We can identify the singlets as spinless fermions, the unpaired impurities as hard core spin 
1/2 bosons hopping on the lattice. Let us consider a system described by the following 
Hamiltonian: 

h = (1/2) ^ P,{U,blkJ,gj)P, (17) 

where the b fields are bosonic, g fields are fermionic, b fields commute with the g fields, and 
the Gutzwiller projector Pg = n^i[^i t .^-p f,t f,. ]. The basis vectors may be represented 
as follows: 
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One can show that the systems described by Hi and h are isomorphic to each other, by 
verifying the following matrix elements: 

<p\Hi\a>=<(3\h\a>, (19) 

where there is the one-to-one correspondence \a >^^ \a > for the basis vectors. The 
Hamiltonian h is equivalent to the following Hamiltonian: 

h = (1/2) Y: {-t,r)PFFlF,,PF, (20) 

where Pp = ULi Ppii), and Ppii) = (1 - F^^Fi^F^^Fn), and L - A^'^ is the number of the F 
fermions on the lattice. 

With the above identification, we have mapped the chiral Kondo lattice model onto the 
chiral Hubbard model with strong repulsion. Therefore, an infinite number of mutually 
commuting invariants can be obtained for the Kondo lattice model. The wavefunctions and 
the thermodynamics of the system may be read from previous results . Any state vectors 
can be written as 

10 >= E H{X},{Y})l[Fy^^YlF},^^Fx,^\P>, (21) 

{X},{Y} i=l j=l 

where \P >= lli=i F--^\0 >, Q = N^ is the number of g fermions, B = -M + {L - Ne)/2 is 
the number of down-spin b bosons, the amplitude $ is antisymmetric in the positions {Y}, 
while symmetric in the positions {X}. The following Jastrow wavefunctions are eigen states 
of the Hamiltonian: 

$({X}, {¥}) = e(2-Vi)("^=S.^-+'"'^S.^^) []d2^Xi - Xj)Y[d{Yi - Yj)lld{Xi - Y,), (22) 

i<j i<j i,j 

with d{n) = sm{iTn/L). The quantum numbers ms,mh are integers or half- integers, which 
make sure of the periodic boundary conditions, satisfying the following constrains 

\mh\ <L/2-{B + Q)/2, 

\mh-m,-L/2\<L/2-{A + Q)/2, (23) 
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with the eigen-energies given by 

E{ms, TTih) = -{2nt/L) [2mh - m, + L/2]Q(l/2). (24) 

The full spectrum of the system takes the following form 

Q Q 

E = -{2Tit/L)[Y,ni + ^ m^](l/2). (25) 

Here, the integers (or half integers ) satisfy the conditions |?t.j| < L/2 — (A + Q)/2, |r7i^| < 
L/2 — (5 + Q)/2, where rii < n^+i and m^ < m^+i. This result shows that the spectrum is 
invariant when changing the sign of t. 

The Jastrow product wavefunctions of the unpaired impurity spins and the singlets are 
typical RVB-type wavefunctions. Various correlation functions of the impurity spins and the 
singlets can be computed exactly, by trivially generalizing Forrester's work to this case. It 
should be remarked that the far away unpaired impurity spins are also strongly correlated 
with each other, because of the fact that only right movers exist in the conduction band. At 
half-filling, the system is obviously an insulator, since each electron forms a singlet at each 
site and the singlets can not hop from one site to another. 

For this chiral Kondo lattice model, the conduction electrons move only in one direction. 
We can anticipate many physical properties for the system even at finite coupling constant 
J . Away from half-filling, one would expect the system to be in a metallic state. At half- 
filling and for sufficiently large J, the system is expected to be insulating. In this model, 
the chirality of the conduction band won't prevent the system from becoming insulating, 
unlike in some other situations, such as in the edges of the Fractional Quantum Hall Effect, 
where the chiral Luttinger liquid won't become localized under any randomness, due to lack 
of backscattering of the quasiparticle |T^. In our case, although the electrons are moving in 



only one direction, the mechanism for localization is very different. The electron always feels 
the exchange interaction of the impurity spin, through the spin exchange interaction. For J 
large enough, each electron will attempt to form a localized singlet with each impurity spin, 
therefore, at half-filling, to transfer one electron from one site to another would break two 
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singlets, causing a charge gap of order of 0{J), and the system would be in an insulating 
state. At J = 0, the system is a simple Fermi liquid. One would thus expect that there 
exists a critical coupling Jc, where the system exhibits metal-insulator phase transition at 
half-filling. An interesting question is whether any infinitesimal small J would drive the 
conducting band to an insulating state at half- filling, i.e. Jc = 0^. If Jc 7^ 0+, one would 
expect Jc ~ |t| by dimensional analysis, and the system is metallic for < J < Jc, while it 
becomes insulating for Jc < J. Metal-insulator phase transition would also occur at half- 
filling when changing \J\ for ferromagnetic interaction between the conduction electrons and 
the impurity spins. Further work is necessary to locate the critical coupling Jc. It might 
also be interesting to see whether the model at finite J belongs to Jastrow-integrable type. 
In summary, we have obtained an infinite number of constants of motion for the one 
dimensional chiral Hubbard model in the strong interaction limit U = 00. We have also 
shown that this model is equivalent to the one dimensional chiral Kondo lattice model at 
J = +00. It seems that the integrability condition might be investigated for finite on-site 
energy, using a similar approach. However, we have not succeeded in doing so for finite on- 
site energy. It has also seemed that the finite J chiral Kondo lattice model is very probably to 
be integrable, as the conduction electrons move in one direction, only exchanging spins with 
the local moments. In the continuum limit, loosely speaking, the many particle scattering 
matrices are very probable to be factorized into two body ones with a continuum relativistic 
electrons, and the situation of scattering matrices of the electrons off one impurity might 



be similar to the S-matrices of the solution of Andrei and al |13|. Further work, either 



numerically or analytically, is necessary for its solvability evidence. 
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